Landau- Ginzburg Perspective of Finite-Temperature Phase Diagrams of a 
Two-Component Fermi-Bose Mixture 
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Abstract 

We consider a mixture of two-component Fermi and (one-component) bose gases under the repulsive Bose-Fermi and attractive 
Fermi-Fermi interaction. We perform a systematic study of the finite-temperature phase diagrams in the chemical potential 
space, identifying, using the Landau- Ginzburg theory, the features generic to the phase diagrams within the validity of our 
model. We apply the theory to explore the physics of correlated BCS pairing among fermions in a tightly confined trap 
surrounded by a large BEC gas. 

PACS: 67.85.-d, 03.75.Ss, 03.75.Mn 



I. INTRODUCTION 

Fermions, constrained by the Pauli exclusion princi- 
ple, behave very differently from bosons, to which such 
a principle does not apply. In a many particle system, 
the former tend to stay away from each other, while, on 
the contrary, the latter tend to be gregarious. In the 
most fundamental level, fermions (leptons and quarks) 
are the building block for all the matter with mass, while 
bosons serve as the mediator for all the fundamental 
forces in nature. In spite of this vast difference, when 
mixed together at temperatures so low that the de Broglie 
wavelength of the particles becomes comparable to or 
even longer than the interparticle spacing, fermions and 
bosons can conspire to create fascinating quantum effects 
at the macroscopic scale that are of fundamental inter- 
est across a broad spectrum of physics, especially within 
the disciplines of condensed matter and nuclear physics. 
As most substances in nature solidify before the tem- 
perature could reach the regime where the macroscopic 
quantum nature of liquid (or gas) can be manifested, the 
combination of liquid isotopes between 3 He and 4 He |l| 
has remained the only laboratory accessible system un- 
til recently when the rapid technological advancement 
in cooling and trapping of neutral atoms has completely 
turned the situation around. Not only has it resulted 
in a dramatic proliferation of such systems, including 6 Li 
- 7 Li 0, i, 6 Li - 23 Na i, 87 Rb - 40 K[|1, 6 Li - 87 Rb 
@, but more significantly, with ultracold atomic gases, 
important parameters, including the interaction between 
particles of same or different species, can be tuned pre- 
cisely @-Qi| , allowing the physics of Fermi-Bose mixture 
to be investigated in a well controlled manner, in regimes 
possibly well beyond the reach by traditional solid state 
systems. 

Recently, by changing the Feshbach detuning across 
a certain critical point at which all the minority atoms 
pair up with the majority ones to form a molecular Bose 
condensate, the group at MIT [ll| has successfully cre- 
ated, from a two-component Fermi mixture with popula- 
tion imbalance, a quantum gas where Bose molecules are 
mixed with spin-polarized (unpaired majority) fermions. 
In addition to confirming an earlier theoretical predic- 
tion of the existence of a transition from full miscibil- 



ity to phase separation [TJI, LL3[ , using the same system, 
they were able to determine the effective dimer-fermion 
scattering length within a reasonable agreement with the 
prediction made more than 50 years ago 14] but never 
verified experimentally, once again demonstrating that 
the ultracold atom system provides an excellent experi- 
mental platform for testing theories. 

Inspired by this work, instead of one Fermi state as in 
a single-component Fermi-Bose mixture in Ref. 



11] 



consider a two-component Fermi-Bose mixture involving 
a hyperfine state |b) of a bosonic atom with mass tub 
and two equally populated hyperfine states: spin up 
I'f) and spin down |i) of a fermionic atom with mass 
niF- The latter Fermi system when equipped with Fes- 
hbach resonance has been the main source of inspiration 
for much recent excitement in the forefront of ultracold 
atomic physics, due chiefly to the vital role it plays in 
the study of crossover from Bose-Einstein condensation 
(BEC) of tightly bound atom pairs to nonlocal Bardeen- 
Cooper-Schrieffer (BCS) atom pairs. Thus, mixtures of 
bosons with such Fermi systems shall be widely accessible 
to experiments. 

The low temperature physics of the two-component 
model under consideration is dominated by s-wave col- 
lisions, which, given that the Pauli exclusion principle 
prohibits s-wave scattering between identical fermions, 
are parameterized with three scattering lengths: aBB, 
apF, and asF, describing, respectively, s-wave scatter- 
ing between two bosons, between two fermions of op- 
posite spins, and between a boson and a fermion of ei- 
ther spin, assuming that scattering a boson off a spin-up 
fermion has the same amplitude as that off a spin-down 
fermion. In principle, clbb must be positive necessitated 
by the Bose stability in the mixture, while aBF and aFF 
can take both positive and negative values as they are 
not subject to similar constraints. Thus, in spite of 
the restriction on obb, such a model can still represent 
drastically different regimes of physics as far as pairings 
and instabilities are concerned. In the present work, we 
aim to extend the physics of phase and phase separa- 
tion at finite temperature from single- to two-component 
systems where fermions have the opportunity to form 
correlated BCS pairs. As a result, we limit our study 
to systems with repulsive (aBF > 0) Bose-Fermi and at- 
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tractive (o,ff < 0) Fermi- Fermi interaction. 

In Sec. UH we review, within the framework of mean- 
field theory, the path integration formulation of the finite- 
temperature thermodynamic potential and derive from it 
the Landau-Ginzburg expansion. In contrast to the ex- 
isting works for the two-component model [Hl-Oj} , which 
have various purposes, we focus on producing phase dia- 
grams in the space made up of chemical potentials, which 
are intensive statistical variables that must remain in- 
variant among the separated phases and hence uniquely 
define a phase separation This is to be contrasted 
to spaces where coordinates are served, for example, by 
particle number densities, of which separated phases in 
a phase separation have different values [H, EH]. In Sec. 
Ill, we apply the analytical intuitions derived from the 
Landau-Ginzburg theory to identify as well as to clar- 
ify the features that are generic to the finite-temperature 
phase diagrams within the validity of our model. The 
utility of such phase diagrams can be most easily appre- 
ciated when one wants to map out the particle number 
densities for a trapped model where the local density ap- 
proximation holds. An example will be shown in Sec. 
IVl which illustrates how the surrounding large BEC af- 
fects the physics of pairing among fermions in a tightly 
confined trap. Finally, we provide a short conclusion 
and discussion in Sec. W\ 



II. MEAN-FIELD THERMODYNAMIC POTEN- 
TIAL AND LANDAU-GINZBURG EXPANSION 

The thermodynamic properties of our model can be de- 
scribed by the partition function Z , which is a functional 
integral over both the complex fields for bosons: ipk, b (t) 
and ip£ B (r) and the Grassman fields for fermions: 

■0k, a ( T ) and V'kjtr (j) in the momentum (Tik) and imag- 
inary time (r) space 20]. In this work, we limit 
our study to the regime of temperature far below Tb 
— 2irh 2 [ub/C (3/2)] 2 ^ 3 /(tub &b), where tib is the Bose 
atom number density, ks the Boltzman constant, and 
C (x) the Riemann-Zeta function. In this limit, bosons 
are virtually all condensed to the zero momentum mode, 
and the standard symmetry breaking ansatz, ipk,B (t) = 
V'o.B + ^k/o.B (t), is therefore applicable, where ipo,B is a 
classical field (not part of the path integration variable) 
for condensed particles, and </>k, b (t) is a bosonic field 
for condensate excitations or simply phonons for a dilute 
gas where the usual Bogoliubov approximation holds [l6[ . 
Further, for pedagogical reason, we opt to first ignore all 
the phonon degrees of freedom, justifying, however, in 
the end of the paper, that the ensuing formalism, when 
appropriately modified, can be applied for a large class 
of Fermi-Bose mixtures where the effects of phonons are 
included. 

The grand partition function Z, at tempera- 
ture T(=l/kB/3), boson chemical potential /is, and 
fermion chemical potential fip, then reads Z = 
J D [tp, if\ exp {-S[n B , ip,ip]/n), where n B = \iPo.b\ 2 /V 



is the boson number density, V the total volume, 5* the 
action given by 

S = h(3V \ ^^-n 2 B - flBTlB 
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In arriving at Eq. ([1}, we have incorporated the inter- 
action between fermions and bosons into Eq. (fTJ) via 
£k = ek — /i in terms of the effective fermion chemical 
potential fj, = [if — 9bf^ B j where ek = (hk) 2 /2itif is 
the kinetic energy of a fermion with a momentum fik, 
and qbb = 47r?i ass/ms and gsF — inh clbf I^bf 
measure the strengths of respective s-wave scatterings 
with ttlb f = 2m biti p / (tub + mj?), and we have also 
modeled the interaction between fermions of opposite 
spins with the parameter U, which, in the actual cal- 
culation, will be replaced in favor of the physical param- 
eter Uq (= 4:Trh 2 aFF /"if) via the standard renormaliza- 



tion relation U 
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In anticipation of the BCS pairing as a result of an at- 
tractive fermion-fermion interaction, we follow the stan- 
dard procedure in which we first introduce the auxil- 
iary bosonic fields A and A* (assuming they are uni- 
form and static) and then apply the Hubbard-Stratonovic 
decomposition to change the partition function into 
Z = fD$,i;,A*,A]exp(-S>[n B ,jj,i;,A*,A}/ti), 
where the action 



S' = hp 
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is expressed in terms of Nambu spinor ipk,iu> n = 
(ipk,tu> n ,t,4>-k,-iu> n ,i) T , with V , k,iw„,t or i being the field 
components in the imaginary (Matsubara) frequency iuj n 
space. 

Finally, by integrating out fermionic fields and carrying 
out a summation over the Matsubara frequency, we ar- 
rive, within the saddle point approximation, at the grand 
thermodynamic potential (density) [SI = — In Zj ((3V)] 
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for a mixture with Bose particle number density ns and 
Fermi particle number density 
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where = y/ ^ + A 2 is the quasiparticle energy, and 
/ (e) = \/{e e / kBT + 1) is the standard Fermi-Dirac dis- 
tribution. Alternatively, we can write Q p defined in Eq. 
(|3]) as a Landau-Ginzburg expansion in terms of the order 
parameter A according to 



Q.f = a Q + a 2 A 2 
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Equations ([3} and ([5]) serve as the foundations for our 
studies below. 
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III. PHASE DIAGRAMS IN CHEMICAL POTEN- 
TIAL SPACE 

In this section, we aim to gain a systematic under- 
standing of the phase diagrams at finite temperature. 
An immediate obstacle to our goal is the existence of 
a large number of free parameters in our model. In or- 
der to alleviate this difficulty, we introduce a parameter 
jimp which has the physical meaning 



system generalized from Ref. [12fl . in which several key 
parameters are scaled according to 



that AnfJ° equals the Fermi energy cp, and adopt a unit 



9bf^ - _ 9 5 bfVb _ _ 9% F ^F 

9BF n B - _ 9% F n F r _ n -~ k 

9% B A^ nF - g% B A^ k - 
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9% B A 6 

n B = 



9bf^ 
9bbA' 



Additionally, A, \i, fc^T, and £k are scaled same as 
fj,p, and the interaction between fermions of opposite 
spins [Uq = Uq/ {g B p / 9bb)~\ is measured relative to 
9bf/9bb - the magnitude of induced fermion-fermion 
interaction due to the density fluctuation of Bose con- 
densate [lij]. As can be seen, the three original free 
parameters (jjbf 9bb ^ Uq) are now reduced into a single 
scaled parameter Uq, making this unit system particu- 
larly suitable for studies aimed at identifying features 
generic to all the phase diagrams. To gain a quali- 
tative understanding of this unit system, we note that 
for a mixture of 87 Rb and 84 Rb with cibb = 100ao and 
a BF = 550a where a is the Bohr radius [2l[, the unit 
is 1.61xl0" 49 J m 3 for Uq, (2.15xlCT 30 J, 3.84xl(T 31 
J) for (fj,p, (j,b), (1.34xl0 19 m^ 3 , 7.49xl0 19 m~ 3 ) for 



FIG. 1: (a) is the illustration designed to aid the understand- 
ing of the zero-temperature phase diagram in (b) where Uo =- 
1.126. Units are defined in the text. 



{rip, ns), and 156 nK for T. (However, for notational 
simplicity, same symbols will be used to stand for the 
scaled variables throughout the rest of the paper.) 

FiguresQJb) and[5]showcase the features that are char- 
acteristic of both the zero- and the finite-temperature 
phase diagrams for a typical system within the validity 
of our model. They are constructed numerically by mini- 
mizing the potential in Eq. ([3]) with respect to and A, 
which amounts to analyzing the saddle point equations 
and the relevant Hessian matrix [l(| [TtJ , as well as com- 
paring, in the coexistence regions, the thermodynamic 
potentials of the possible phases listed below: 

(V) for an empty phase where both ns = and np = 

0, 

(B) for a pure boson phase where ns ^ and np = 0, 
(N) for a normal state of Fermi gas where ns = 0, 

A = 0, and np ^ 0, 

(S) for a superfluid Fermi gas where iig = 0, A ^ 0, 

and uf ^ 0, 

(BN) for a normal Fermi-Bose mixture where ns i= 0, 
A = 0, and np ^ 0, and finally 

(BS) for a superfluid Fermi-Bose mixture where ^ 
0, A ^ 0, and n F ^ 0. 

In order to capture the physics, making the study of 
phase diagrams more illuminating, we complement our 
numerical study with an analysis founded on the Landau- 
Ginzburg theory in the low temperature limit T « fi. 
In this limit, [Eq. (|8c[) below] is always positive, and 
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we only need to consider the Landau-Ginzburg expansion 
up to the fourth order in A, 
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where the integrals for the Landau coefficients in Eqs. 
© can be evaluated explicitly with the results [22| 



ao 



a2 



1 



a! 

4 V/' 



2 
5 

3 1/2, ^ 

4^ ln 8^' 



„ 21 , ( 1/2C(3) 



(8a) 

(8b) 
(8c) 



with 7 0.577 being the Euler's constant. Further, in 
the same limit, we can ignore the thermal population in 
comparison with the population inside the Fermi sphere 
so that the Fermi density can be approximated as 



np = 



on 
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TT 2 T 
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Equations JSj) and ([9]) form the backbone of the analytical 
discussions that we carry out below. 

Consider first the mixed phase where both ns and Uf 
have finite values. We begin with the saddle-point equa- 
tion for ub 



nB - Mb + M 



3/2 



0, 



(10) 



which, when the use of [i = /if — ?*b is made, leads to a 
cubic equation foryTI 



(VJ 1 ) 3 ~ (Vtf + (mf - Mb) = 0. 



(11) 



The positive root to Eq. (jlll) that lies below 2/3 (oiyfjl 
< 2/3) represents a local stable mixed phase, where the 
inequality is derived by subjecting Eq. (0 to the local 
stability criteria that d 2 il/dn B be positive. A simple 
analysis then shows that Eq. (|TT|) supports the local sta- 
ble mixed phase in the chemical potential space between 
the line [labeled as 1 in Fig. [Ha)]: 



Mb = MF, 



(12) 



and the parallel line [labeled as 2 in Fig. QJa)] shifted 
with respect to Eq. {12} by 4/27: 



fi B = Mf - 4/27. 



(13) 



Next, consider the pure Bose phase described by the 
bosonic potential Ub — n B /2 — fiBns and the pure Fermi 

9 5 /2 

phase described by the fcrmionic potential £lp — — f Mf ■ 
In the region that supports the pure Fermi (Bose) phase, 
the effective chemical potential for fermions fip ~ n B 



(bosons fiB — np ) must be negative. As a result, the 
pure Bose phase exists above the line in Eq. (fl2")) . and 
the pure Fermi phase exists below the curve [labeled as 
3 in Fig. [Ha)] 



Ms 



3/2 

Hf ■ 



(14) 



It is then clear that the pure Fermi phase coexists with 
the mixed phase in the region below curve 3 but between 
lines 1 and 2, and it also overlaps with the pure Bose 
phase in the upper right region below curve 3 and above 
line 1. 

So far, nothing has been said regarding the nature of 
the Fermi gas component. To answer this question, we 
go to the saddle point equation for A 2 



A 



«2 (/yT) 

a 4 (fi,T) ' 



(15) 



Evidently, because «4 (/i, T) is always positive, only when 
a-i (/i, T) < does BCS order or superfluidity occur. Let 
Hp be the root to the threshold condition 



Q 2 



(4> T 




0. (16) 



In the limit where T goes to zero, fip also goes to zero, 
and hence as expected under the attractive Fermi-Fermi 
interaction, fermions at T = either exist as a pure S 
state or mix with bosons to form a mixed BS state. 

Finally, by comparing the energies of the phases in each 
coexistent region, we change Fig. [Ha) i n ^ the phase 
diagram in Fig. [ljb), which features a tricritical point 
(labeled as *) at which the transition between BS and 
S changes from the second- to first-order type, and a 
critical point (labeled as x) at which the second-order 
transition line between B and BS and the two first-order 
transition lines between S and, respectively, BS and B 
meet. The explanations are provided as follows. The 
first-order transition line (dotted) between B and S is 
given by 



Mb = \/4/Ei/J 4 ', 



(17) 



which is obtained by setting Qb = ^f- The interception 
between Eq. (fT2|) and Eq. (fl7|) defines the critical point 
(fi F = n B = 25/16). The comparison between fi and 
determines the first-order transition between BS and 
S, which starts from the critical point x and ends at the 
tricritical point * with (fip = 4/9, fi* B = 8/27) at which 
jib = dfl/dris — d 2 fl/dn B = 0. 

Note that for a single-component model with the repul- 
sive Bose-Fermi interaction, the phase diagram at T = 
has the same structure as FigOJb) except that S and 
BS are replaced, respectively, with N and BN, owing to 
the absence of pairing mechanism for fermions 

mi. As 

stressed in the introduction, an advantage of the phase di- 
agram in the chemical potential space is that each phase 
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FIG. 2: The phase diagram (a) at T = 0.001 and (b) at 
T — 0.25. £/o is same as in Fig. [T] and units are defined in 
the text. 



separation manifests itself as a unique coexistence curve, 
making its identification particularly transparent. For 
example, Fig. [1] (b) automatically rules out phase sep- 
arations other than those between a pure Fermi and a 
pure Bose phase and between a pure Fermi and a mixed 
phase. Thus, there is no need to hypothesize and analyze 
the existence of, for example, a phase separation between 
a pure bosonic and a mixed phase or that between two 
different mixed phases, as did in works where phase anal- 
ysis was performed in space made up of particle number 
densities [II GJ]. 

At (low) finite temperature [Fig. E[a)], the pure Fermi 
phase is divided into N and S by the vertical (second- 
order transition) line [Lf — ftp, where fip increases with 
temperature according to Eq. (1161) , and the mixed phase 
is split into BN and BS by the (second-order transition) 
line, [lb — [If + tip — fip, which begins at point # and 
ends at point x , with rip being the Fermi number density 

[Eq. (HJ] at A = and fi = fif,. As the temperature 
increases, both the N and BN region expand while S and 
BS contracts. A further increase in the temperature so 
that up" lies beyond fip at the tricritical point results in a 
complete elimination of the phase space for the superfluid 
Fermi-Bose (BS) mixture as shown in Fig. [2fb), where 
the normal state Fermi-Bose mixture remains as the only 
mixed phase. 

At this point, we recall that our analysis rests upon an 
assumption that the term 7r 2 T 2 /8/Lt 2 in Eq. © be much 



less than 1, which is therefore applicable only to situa- 
tions where the Fermi surface is well defined (relative to 
the temperature of interest). At zero temperature, the 
Fermi surface is fixed by /j = 4/9 for line 2 (and the tri- 
critical point), but is not defined for line 1 where \i = 0; 
line 1 holds only for T — 0. As a result, the line dividing 
BN and BS in Fig. [21(a) (which corresponds to line 1 at 
zero temperature in Fig. [T]) is far more sensitive to the 
temperature increase than the first-order transition line 
(and the tricritical point). 

The line between B and BN requires a special explana- 
tion. At zero temperature, due to the Fermi-Dirac distri- 
bution, the chemical potential for a normal state of Fermi 
gas must be positive, and under such a circumstance, the 
pure Bose phase (B) is well-defined - the space above 
line 1, where the chemical potential is negative /i < 0, 
is defined as the phase B. As temperature increases, 
the Fermi surface becomes less sharply defined, and so 
does the distinction between B and BN. In fact, at finite 
temperature there is nothing to prevent the chemical po- 
tential for fermions from becoming negative, just as for 
distinguishable particles that follow the Maxwell's distri- 
bution. In our study here, the line between B and BN is 
defined in such a way that the Fermi density along this 
line is 10~ 12 , a small threshold artificially introduced to 
indicate that the region in B is virtually free of fermions. 
But, we stress that in contrast to the line between BN 
and BS across which there is a second-order phase tran- 
sition due to the emergence of the order parameter A 
which breaks the symmetry, there is no real phase tran- 
sition for fermions across the line between B and BN. 



IV. AN APPLICATION: TIGHTLY TRAPPED 
FERMIONS EMBEDDED IN A LARGE BEC 

In practice, phase diagrams come in different forms, 
depending on the application at hands. A complete 
mapping of the phase diagram including all the criti- 
cal points in the chemical potential space as we have 
done in the previous section is significant, not only be- 
cause it provides a roadmap for constructing phase di- 
agrams in other spaces, but also because it can greatly 
facilitate both the identification and the interpretation 
of the important features emerging from these diagrams. 
As an application, we consider, in this section, a double 
species experiment in which spherically symmetric har- 
monic trap potentials, Vb,f (t), are tuned to differ in 
such a manner that fermions are tightly confined inside 
a much larger BEC, and use it as a model to explore 
how the surrounding bosons affect the physics of Fermi 
pairing. This model is the two-component analog of the 
single-component system experimentally realized by the 
MIT group in 2002 [4]. It also represents an example 
in which the traps are oppositely arranged compared to 
those adopted to implement the idea of using a small 
trapped BEC to probe the properties of a large Fermi 
gas component [T3, HH . 
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In our study below, we assume that the Fermi com- 
ponent, although small compared to the BEC, is still 
sufficiently large so that the local-density approximation 
is directly applicable. Thus, for fermions, we introduce 
a local chemical potential pp (f) — fJ-F — Vf if), in ad- 
dition to a global chemical potential pp, which is fixed 
by the total particle number. In contrast, for bosons, 
owing to its large size, to a good approximation, we can 
regard its chemical potential, within the distance scale in 
the order of the size of the Fermi gas, to be a constant, 
Pb, independent of the variation of the radial distance 
r. 

Figure [3] (a) displays the phase diagram for a homoge- 
neous system in the T — pp space for ps = 0.35 [slightly 
smaller than p* B in Fig. Eta)]. At a given temperature 
and within the local-density approximation, a potential 
profile Vf (r) along the radially outward direction in the 
real space induces an image pp (r) in the form of a ver- 
tical line moving down an arrow as indicated in Fig. [3] 
Then, by traversing from the tail to the head of the ar- 
row and translating each point [/if (?")] 011 the vertical 
line into the density via the mapping p F (fip (r)), we can 
construct, from the center to the edge of the trap, the 
density distribution pp (r) = p F (pp (r)), where p F (x) 
calculates the Fermi density for a homogeneous phase at 
the chemical potential x. Using the rules just outlined 
and given the locations and the sizes of the vertical lines 
indicated in Fig. [2t a ), we can easily arrive at the follow- 
ing qualitative conclusions. 

(a) At sufficiently high temperature (the arrow on the 
right), the density profile consists of a normal phase (N) 
core surrounded by a (spherical) shell of a BN mixture. 

(b) As temperature is lowered (the arrow in the mid- 
dle), superfluidity (S) emerges from the core with N being 
sandwiched between S and BN. 

(c) As temperature is further reduced (the arrow on 
the left), N disappears and substituted in its place is the 
Bose-superfluidity mixture (BS). 

(d) A further reduction in the temperature completely 
eliminates BN, creating a density profile, reminiscent of 
the high temperature case in (a), except that N and BN 
in (a) are now replaced with S and BS, respectively. 

An intriguing aspect of such a mixture is that both 
critical temperatures and Fermi profiles depend strongly 
on hb - the chemical potential of the surrounding BEC. 
As pb can be independently tuned in experiments, we 
have at our disposal a convenient tool to selectively ac- 
cess the parts of phase diagrams, that are physically in- 
teresting. Suppose that we want to access the first-order 
phase transition. We can do so simply by increasing ps 
beyond a certain threshold value. An example is given 
in Fig. [2tb) where ps is fixed to 0.45, a value slightly 
higher than p* B in Fig. Et a )- As can be seen, part of the 
transition line dividing pure and mixed states changes its 
nature from second order in Fig. [3ta) to first order in Fig. 
[3tb). Thus, although the system can undergo a similar 
set of phase transitions as in Fig. [3t a ) [at the expense 
of a higher Fermi chemical potential (or density) and a 




FIG. 3: The phase diagram in the T — pf space when the 
chemical potential for the bosons is fixed to (a) ps = 0.35 
and (b) ps = 0.45. The features in (a) and (b) can be easily 
deduced from the phase diagrams in the previous section. Uo 
is same as in Fig. [1] and units are defined in the text. 



smaller BS region], a density discontinuity is expected to 
emerge from Fermi profiles on a sphere that divides the 
two phases sharing the same first-order transition line. 

V. CONCLUSION 

A unique advantage of a cold atom system is that 
its system parameters can be precisely tuned, allowing 
it to access a larger regime of phase diagrams. Fur- 
ther, availability of detection techniques, such as ab- 
sorption laser imaging of densities and radio-frequency 
(RF) spectroscopy [2J], makes the experimental deter- 
mination of such phase diagrams in fine detail possible. 
In this work, we have performed a systematic study of 
the finite-temperature phase diagram in the chemical po- 
tential space for a two-component Fermi-Bose mixture 
with attractive Fermi- Fermi and repulsive Fermi-Bose in- 
teraction. Using a combination of scaling and Landau- 
Ginzburg theory, we have identified, within the frame- 
work of mean-field theory, a set of features generic to the 
phase diagrams for such mixtures. Further, we have ap- 
plied the theory to explore the physics of pairing among 
fermions in a tightly confined trap surrounded by a large 
BEC gas. 

Finally, we comment that so far, our analysis is 
based on the thermodynamic potential in Eq. ([3]) 
derived when all the phonon fields, 0k. b (t), are ig- 



G 



nored. In more realistic situations, we need to re- 
tain the phonon fields, 4>k,B (t), which, after being ex- 
plicitly integrated out under the Bogoliubov approxi- 
mation, is shown to induce an attractive interaction, 

U md (q) = -9 2 bf /9bb/ 



1 + 



be- 



[hq/ \/AmBgBBn B ) ' 
tween fermions when the retardation effect is ignored 
[l|| . The net effect of </>k,s (t) is then to change 
Uq to U e M = Uq — 9%f/9bb, for the s-wave scattering, 
which is expected to dominate all the other partial- wave 
scatterings. Thus, the phase-diagram features presented 
in this work shall remain qualitatively true not only for 
two-component Fcrmi-Bose mixtures with direct Fermi- 
Fermi attraction (Uq < 0), but also for those with di- 



rect Fermi-Fermi repulsion (Uo > 0) but attractive effec- 
tive Fermi- Fermi interaction (U e ^ < 0) [26]; the latter 
case turns out to cover several important systems of cur- 
rent experimental interest, including 40 K- 87 Rb and 6 Li- 
7 Li M- 
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